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The unconventional critical behavior near magnetic quantum phase transitions in various heavy- 
fermion metals, apparently inconsistent with the standard spin-density-wave scenario, has triggered 
proposals on the breakdown of the Kondo effect at the critical point. Here we investigate, within 
one specific scenario, the fate of such a zero-temperature transition upon coupling of the electronic 
to lattice degrees of freedom. We study a Kondo-Heisenberg model with volume-dependent Kondo 
coupUng - this model displays both Kondo volume collapse and Kondo-breakdown transitions, as 
well as Lifshitz transitions associated with a change of the Fermi-surface topology. Within a large- A" 
treatment, we find that the Lifshitz transition tends to merge with the Kondo volume collapse and 
hence becomes first order, whereas the Kondo breakdown transition remains of second order except 
for very soft lattices. Interesting physics emerges at the zero- temperature endpoint of the Kondo 
volume collapse: For electrons in two space dimensions, this endpoint is located at the Lifshitz line 
for a large range of parameters, thus two continuous quantum phase transitions coincide without fine 
tuning. We analyze the effective Landau theory for such a situation and discuss critical exponents. 
Finally, we relate our findings to current heavy-fermion experiments. 



I. INTRODUCTION 

The past decade has seen growing interest in continu- 
ous quantum phase transitions. The associated finite- 
temperature quantum critical regime is characterized 
by unconventional dynamics and thermodynamics; in 
particular, quantum criticality can profoundly modify 
the metallic properties of both d- and /-electron ma- 
terials. Heavy-electron metals have played a particu- 
larly important role in the study of quantum critical- 
ity, for these materials lie at the brink of antiferromag- 
netic (AF) instabilities, and they are readily tuned to 
an AF quantum critical point (QCP). The occurrence 
of this QCP is commonly attributed to a competition 
between the Kondo screening of the local moments and 
a tendency to AF ordering due to an inter-moment ex- 
change interaction, as first proposed by Doniach.^ In- 
terestingly, various properties near the critical points 
of materials like CeCue-^Au^; and YbRh2Si2 appear 
to be inconsistent with the standard spin-density-wave 
(or Landau-Ginzburg- Wilson) scenario of the magnetic 
phase transition.^ Most prominent here are the singular 
behavior of the specific- heat coefficient and £^/T-scaling 
in inelastic neutron scattering experiments.^ These have 
led to suggestions that the critical point is not simply a 
magnetic instability of fermionic quasiparticles, but that 
in addition the Kondo effect, being responsible for the 
formation of these quasiparticles, breaks down.^^ 

Assuming that the Kondo breakdown is the primary 
source of the critical singularities, Senthil et al.^'^ have 
suggested a model where the Kondo breakdown occurs 
as a topological quantum phase transition between two 
paramagnetic phases, namely a heavy Fermi liquid (FL) 
and a so-called fractionalized Fermi liquid (FL*), and 
ordered magnetism is regarded as an epiphenomenon. 
As in the Doniach picture, the formal control parame- 
ter of the transition is the ratio Tk / Jh between Kondo 



temperature and Heisenberg exchange of local moments. 
It is, however, assumed that the local moments form 
a spin-liquid state without broken symmetries (instead 
of the standard local-moment antiferromagnetism) in 
the regime Tk <C Jh- Such a spin liquid may result 
from strong quantum fluctuation, induced e.g. by frus- 
trated RKKY or exchange interactions between the local 
moments.^ 

As phase transitions in heavy-fermion metals are often 
tuned by external or chemical pressure, the coupling of 
electronic to lattice degrees of freedom can play a non- 
trivial role for the overall shape of the phase diagram. 
This aspect motivates our paper. A well known but 
spectacular example are the pressure-induced volume- 
collapse transitions in the trivalent rare earth metals,^ 
like the transition between a and 7 Ce. In this ex- 
ample, a line of first-order transitions is found in the 
pressure-temperature plane where the unit cell reduces 
its volume by about 15%. An order of magnitude change 
of the Kondo temperature accompanies this transition. 
These aspects have led to the Kondo volume-collapse 
model, ^^'^^ which explains the transition based on a vol- 
ume dependence of the Kondo exchange coupling. 

In this paper, we shall study a model which contains 
the ingredients both for the Kondo-breakdown quantum 
phase transition as well as for Kondo volume-collapse 
physics, with the goal of obtaining a unified phase di- 
agram of these phenomena. As both transitions can in 
principle be tuned by pressure, with the inter-atomic 
distances strongly influencing the Kondo energy scales, 
an intricate interplay of electronic and lattice degrees of 
freedom may be expected. Of particular interest will be 
whether the Kondo-breakdown transition of the type pro- 
posed in Ref. 7 survives as a second-order quantum tran- 
sition once a coupling to the lattice is taken into account. 
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A. Summary of results 

Our microscopic analysis is based on a fermionic large- 
N treatment of a Kondo-Heisenberg model in both two 
and three spatial dimensions, with the Kondo exchange 
being linearly coupled to static lattice distortions. In the 
scenario of Ref. 7, the Kondo-breakdown transition be- 
tween the FL and U(l) FL* phases occurs in the presence 
of two Fermi surface sheets; as a result, a Lifshitz tran- 
sition generically occurs upon tuning the ratio Tk/Jh 
within the heavy Fermi-liquid phase. At such a Lifshitz 
transition, well known from weakly interacting Fermi 
systems, a topological splitting of the Fermi surface oc- 
curs, here leading to a transition between Fermi liquids 
with one and two Fermi sheets (dubbed FLi and FL2, 
respectively) . 

The physical behavior of our system depends on both 
the ratio Tk/Jh and the electron-lattice coupling (or 
the lattice stiffness). Zero-temperature phase diagrams, 
as function of pressure (which linearly couples to Tk) 
and Ji^, are shown below in Figs. 2, 3. Interestingly, the 
above-mentioned Lifshitz transition is easily driven first 
order upon inclusion of lattice degrees of freedom, and is 
then accompanied by a jump in lattice volume and Kondo 
temperature. In contrast, the Kondo breakdown transi- 
tion FL2-FL* remains of second order except for very soft 
lattices - in the latter case, the Fermi-liquid breakdown 
occurs from FLi to FL* as a strong direct first-order tran- 
sition. For small J//, the behavior is dominated by the 
interplay of Kondo and lattice physics, and soft lattices 
admit a separate Kondo volume collapse transition (as in 
earlier work) . This first-order transition merges with the 
Lifshitz transition at larger Jh- 

For a large range of lattice stiffness values, a zero- 
temperature endpoint of the Kondo volume collapse oc- 
curs in the p-Jn phase diagram. Remarkably, in two spa- 
tial dimensions this endpoint is quite generically located 
on the Lifshitz transition line without further fine tuning. 
This coincidence of two continuous transitions can be fur- 
ther analyzed using an effective theory for coupled Lif- 
shitz and lattice order parameters. It turns out that the 
electron-lattice coupling leads to non-perturbative effects 
which are ultimately responsible for the "non-Landau" 
transition behavior. 



B. Relation to earlier work 

Our work merges ideas on the Kondo volume collapse 
and those of quantum critical behavior near magnetic 
phase transitions in heavy- fermion compounds. Further- 
more, earlier works on quantum Lifshitz transitions play 
a role as well. 

Historically, two main theoretical approaches have 
been advanced to describe the 7 ^ a transition of el- 
emental Cerium, namely the Mott transition scenario-"^^ 
and the Kondo volume collapse model. ^^'^^ Common to 
both scenarios is that the / electrons are more localized in 



the 7 phase. More recently, the modern dynamical mean- 
field theory (DMFT) in combination with realistic band 
structure calculations (LDA+DMFT) was employed, and 
reasonable agreement with experimental data has been 
achieved. The relation between the Kondo volume col- 
lapse and quantum criticality was only discussed recently 
in Ref. 15, where the situation of the endpoint suppressed 
to T = was considered. This work did not include com- 
peting exchange interactions between local moments, and 
hence did not contain the physics of magnetic or Kondo 
criticality. 

The breakdown of the lattice Kondo effect at the QCP 
in heavy fermions was proposed in Refs. 4,5. In par- 
ticular, the approach of Si et al.^ used an extension of 
DMFT to obtain a scenario of "local quantum critical- 
ity", where the onset of magnetism in two dimensions 
suppresses Kondo screening right at criticality. A some- 
what different viewpoint was put forward in Refs. 6,7, 
where Kondo criticality was discussed in the absence of 
ordered magnetism. This required the existence of a 
novel phase, a fractionalized Fermi liquid (FL*), which is 
a metallic spin-liquid phase with emergent gauge struc- 
ture. Central to the Kondo-breakdown transition^'^ is 
that the / electron subsystem becomes insulating, hence 
the transition can be described as an orbital-selective 
Mott transition. ^^'^^ Ref. 7 studied the mean- field theory 
to be used below, but the coupling to lattice degrees of 
freedom was not discussed, and the Lifshitz transition in 
the FL phase (which will turn out to be important in our 
analysis) was mentioned, but not investigated in detail. 
Similar mean- field theories also appear in Refs. 18,19, but 
not with focus on the Kondo breakdown transition. 

Lifshitz transitions in correlated electron systems, with 
topological changes of the Fermi surface, have been dis- 
cussed at length in Ref. 20. However, this analysis did 
not contain the coupling of the electrons to other collec- 
tive degrees of freedom. 



C. Outline 

The body of this paper is organized as follows: In 
Sec. II we formulate the Kondo-Heisenberg model supple- 
mented by electron-lattice coupling and discuss its large- 
TV solution and the resulting zero-temperature phases. 
Fully self-consistent numerical calculations of the zero- 
temperature phase diagrams in the parameter regime of 
pressure and inter-moment exchange are presented after- 
wards. Particular attention is paid on the role of di- 
mensionality and geometry of the underlying lattice. We 
discuss the existence of discontinuous volume collapse 
transitions that occur in a wide range of our numeri- 
cally obtained phase diagrams. In Sec. HI we focus on 
the interplay of Lifshitz and volume-collapse transitions 
in the framework on an effective Landau theory. This 
analysis is in agreement with the mean-field results for 
the Kondo lattice, and furthermore shows that, for two- 
dimensional (2d) electrons, the quantum critical endpoint 
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of the first-order volume collapse transition tends to co- 
incide with the second-order Lifshitz transition. Fluctu- 
ations beyond mean-field theory are discussed in Sec. IV, 
and a brief summary and outlook concludes the paper. 



II. MICROSCOPIC KONDO-LATTICE THEORY 

The starting point of our analysis is the standard An- 
derson lattice model, describing conduction electrons on 
a lattice which hybridize with local atomic /-orbitals on 
the lattice sites. 

ka i(J i 

+ ;4E(^^«"'''''4./i. + ^^-c.) (1) 

^ kia 

The first term describes conduction electrons with band 
filling nc, moving on a two-dimensional (2d) or three- 
dimensional (3d) regular lattice (the former case applies 
to layered systems with weak electronic inter-layer cou- 
pling), £j {U) is the bare / electron energy (Coulomb 
repulsion), Vk the hybridization matrix element, and Af 
the number of unit cells. 

Let us now discuss the influence of volume changes 
induced by external pressure. Experimental data on 
Kondo volume collapse transitions in Ce shows that, 
throughout the transition, the / level occupation and 
the / electron levels show only modest changes. In con- 
trast, the Kondo temperature changes rapidly with vol- 
ume by orders of magnitude. These results motivate us to 
neglect any volume dependence of the Anderson lattice 
model other than that of the hybridization matrix ele- 
ments. (Note that the Kondo temperature depends expo- 
nentially on both the hybridization and the conduction- 
electron bandwidth. Experimentally, both will be influ- 
enced by volume changes. For our qualitative considera- 
tions, taking into account the volume dependence of one 
is sufficient.) 

For most of the paper, we will treat the lattice degrees 
of freedom classically, parameterized by the strain ten- 
sor e. At a transition between two iso-structural states 
a suitable order parameter is the trace e of the strain 
tensor, 
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where Vq is a reference volume. In analogy to Ref. 15, the 
hybridization matrix elements are allowed to depend on 
the local strain by a linearized change of the overlap 
integral Vk = {k\V\f) between a Bloch electron and a 
local / electron with its atomic potential V. Assuming 
a local hybridization, Vk = V ^ we account for the lattice 
coupling be replacing Ti ^ Ti-\- Tic with 
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where 7 is the coefficient of the assumed linear lo- 
cal strain dependence of the hybridization. The last 
term captures the deformation energy, with B being 
the bulk modulus. For instance, in a cubic system, 
B = \{c\i + 20^2), with bare elastic constants and 
^0 21 

Because of the exponential variation of Tk with J^, 
moderate variations of the Kondo coupling Jk with vol- 
ume cause a strong volume dependence of the Kondo 
temperature. This, in turn, yields a nonlinear equation 
of state j9(y) = -ff It = -^I^It- It has been shown^^ 
that models of this type yietd nonlinear p-V isotherms 
similar to the van der Waals theory of the liquid-gas tran- 
sition. The phase diagram of the model from Eq. (1) and 
(3) has been discussed in Ref. 15, and a zero-temperature 
volume collapse within the heavy- Fermi liquid phase has 
been found to occur below a critical value B^ of the bulk 
modulus. 

In this paper, we shall discuss the interplay of volume 
collapse transitions with a possible breakdown of Kondo 
screening due to competing inter-moment exchange. To 
simplify the following approximate treatment, we shall 
supplement our model by an explicit Heisenberg-type ex- 
change interaction between the local moments,^' and we 
shall pass from the Anderson to the Kondo model. The 



latter step, justified in the Kondo limit 
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neglects possible charge fluctuations which we assume to 
be of minor influence for the structure of phase diagram. 
Note that we only fix / valence (which is known to not 
vary drastically across a volume collapse transition) to 
unity, whereas we still allow for a transition from "itiner- 
ant" to "localized" / electrons.^ We arrive at a Kondo- 
Heisenberg model of the form 

WkHM = ^^kcl^C^^ + JK{l+ief^Sr'Sr 

k r 



-\-Jh ^ ^ Sr • Sr' 



(4) 
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where Sr = ^cl^acra'Cra' local spin density of the 

conduction electrons. We neglect the influence of lat- 
tice distortions on inter-orbital magnetic exchange, for 
two reasons: (i) The exponential dependence of Kondo 
temperature on the strain likely dominates over other in- 
stances of electron-lattice couplings, (ii) The physics is 
primarily influenced by the ratio Tk/ Jh^ and for quali- 
tative considerations the strain dependence of one of the 
parameters is sufficient. 

In the following, the zero-temperature phase diagram 
of the model Wkhm (4) will be determined within a 
fermionic mean-field approach,^ which allows to capture 
both the volume collapse and the Kondo breakdown tran- 
sitions. An implicit assumption is that quantum fluctu- 
ations of the local moments are large, such that a spin 
liquid replaces ordered magnetism for Tk <C Jh • Techni- 
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cally, this assumption determines the choice of the mean- 
field decouphng of Jh to be described below. 



A. Fermionic large- A^^ theory 

The model (4) cannot be solved analytically as it 
stands. Instead, much progress can be made by solving it 
in the large- approach which has been extensively used 
in the study of the Kondo lattice and also allows to deal 
with the magnetic exchange term in Eq. (4).^^'-'^^ 

Using a fermionic representation of the local moments 
allows us to capture phases both with and without Kondo 
screening. Importantly, for dominant exchange inter- 
action Jh^ one obtains a spin-liquid or valence-bond 
ground-state,^'^ which replaces the metallic antiferro- 
magnet conventionally found in this regime.^ Here, we 
will restrict our attention to phases without broken lat- 
tice symmetries, hence tuning the ratio Jh/Tk drives 
a transition between a heavy Fermi liquid (FL) and a 
fractionalized Fermi liquid (FL*) - in the latter phase 
the conduction-electron and local-moment subsystems 
are decoupled, and the / electrons form a Mott-insulating 
spin liquid state. The quantum phase transition between 
FL and FL* is accompanied by a jump in the Fermi 
volume.'' 

The large- A/" approach to the Kondo-Heisenberg model, 
Eq. (4), extends the local- moment symmetry to the 
SU(A/') symmetry group. The N'^ — 1 components 

(5'J^/)i<o-,o-'<iv of a local moment S{r) are represented 
in terms of neutral fermions fra with a local constraint: 
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The interaction terms of Eq. (4) now become quartic. 
These are decoupled introducing Hubbard- St r at onovich 
fields br{r) conjugate to J2a fla^ra and Xij(r) conjugate 
to fiafja and constraint is enforced via Lagrange 
multipliers Ar(r). The action corresponding to the rep- 
resentation (5) of Eq. (4) becomes extensive in N by 
rescaling the coupling constants as Jk ^ Jh ^ 
In the N ^ oc limit the physics is controlled by a sad- 
dle point at which the Bose fields condense. We assume 
a spatially homogeneous mean-field state^^'^^ where br^ 
Ar-, and Xij take constant values bo, A, and xo- The 
mean-field Hamiltonian takes the form^^ 
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where we have dropped additional constants. The saddle 
point equations can be written in the following compact 
form: 
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where £kf = A + 2xo X]f*=n n dispersion of the 

/-fermions and Gc, Gf and Gfc are the full conduction- 
electron, /-electron and mixed Green functions, respec- 
tively, obtained from the Matsubara Green functions: 

GciiuJn, k) = GliiuJn f^c- blG^f{iu;n, k), k) 
Gf(iuJn, k) = G^f(iuJn - blG^^(iuJn, k), k) 

boG^iUn.k)G^,{iUn.k) 



Gfc{iuJn,k) 



blG){iuJn,k)G^,{iuJn,k) 



(8) 



Here G^(zcc;n, k) = {iuOn — £k) ^ and G'j{iuJn^ k) = {iuOn — 
Skf)~^ are the non-interacting conduction and / electron 
Green functions. 

The equilibrium lattice strain e minimizes the free- 
enthalpy function 

G(e) =F(e)+pFo6 (9) 
what leads to the additional self- consistency condition 
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which completes the set of mean-field equations together 
withEq. (8). 



B. Zero-temperature phases 

In Eq. (6), the eff'ect of the strain e is just to 
renormalize chemical potentials and the Kondo coupling 
Jk' Thus, there are two qualitatively different zero- 
temperature phases of our mean-field model which where 
already derived in Ref. 7. On the one hand, there is 
the usual Fermi-liquid (FL) phase when ^o, Xo? A are 
all nonzero. (Note that 6o 7^ implies xo 7^ 0-) One 
the other hand, the Kondo hybridization may be zero, 
6o = 0, but Xo 7^ (and A = 0) - this is the FL* phase. 
In this mean-field state the conduction electrons are de- 
coupled from the local moments and form a smalf Fermi 
surface. The local-moment system is described as a U(l) 
spin liquid with a Fermi surface of neutral spinous. We 
expect that xo ^ Jh- At finite temperatures the phase 
transition becomes a crossover - this feature is not cap- 
tured by mean-field theory but requires to consider the 
coupling to the compact U(l) gauge field.'' (Spin liquid 
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FIG. 1: Fermi surface evolution from FL* to FL, where shaded 
areas correspond to occupied states. Left: FL*, with one 
spinon (dark) and one conduction electron (light) sheet. Note 
that the spinon band is hole-like. Middle: FL2, with two 
sheets, where the outer one represents heavy quasiparticles 
with primarily / character. Right: FLi, with one heavy- 
electron sheet. FL2 and FLi are separated by a Lifshitz tran- 
sition where the outer Fermi sheet disappears at critical value 
of the ratio Jh/Jk- The conduction band has a filling of 
Tie — 0.8. (The corresponding band structures are also shown 
in Fig. 4 below.) 



states where the gauge group is broken down to Z2 are 
also possible: Those lead to superconductivity masking 
the FL-FL* transition^ and shall not be considered here.) 

We note that the large- mean-field approach can be 
expected to be qualitatively correct at zero temperature: 
The FL phase is stable w.r.t. fluctuations, as is FL*, pro- 
vided that the U(l) gauge field is in a deconfined phase. 
For a further discussion of fluctuation eflFects we refer the 
reader to Ref. 7. 

The mean-field Hamiltonian can be interpreted as a 
model of non-interacting quasiparticles with dispersions 
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Now consider the FL phase near the transition (small 
&o)- In this case, the two bands derive from the c elec- 
trons (with weak / character) and the / particles (with 
weak c character). For small 60, both bands together 
therefore intersect the Fermi energy at least twice so that 
the Fermi surface consists of two or more sheets (see Fig. 
1). Upon increasing 60, a topological transition is possi- 
ble beyond which the Fermi surface intersects the lower 
band only once and the upper band E^ is empty 
(Fig. 4 below). This topological splitting of the Fermi 
surface is a Lifshitz transition. 

Let us point out that the existence of this Lifshitz tran- 
sition in the FL phase near the Kondo breakdown transi- 
tion can be seen as a consequence of the fermionic nature 
of the spinous in the FL* phase, which then implies the 
existence of two fermionic bands near the FL-FL* transi- 
tion (also beyond mean-field). We will comment on other 
scenarios toward the end of the paper. 
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FIG. 2: Large- ground state phase diagrams of the Kondo- 
Heisenberg model (6) on a 2d square lattice, as function of 
pressure p and inter- moment exchange Jh- The panels a) to 
d) are for bulk moduh 5=0.005, 0.0055, 0.007, and 0.015; the 
other parameters are ric — 0.8, Jk — 1.5 and 7 — 0.05. 
Thin (thick) lines are second (first) order transitions. The 
Lifshitz transition line separates FLi and FL2, the Kondo 
breakdown transition separates FL from FL*, for details see 
text. (The calculations have been at a low T = 0.005, the 
dashed lines are extrapolations of the phase boundaries ob- 
tained from runs at lower T.) 



C. Numerical results: Phase diagrams 

We have solved the mean-field equations (8) and (10) 
numerically for different lattice geometries, namely for 
electrons with nearest-neighbor hopping t on 2d square 
and triangular lattices and for 3d cubic lattices. ^^'^^ The 
geometric frustration (here of the triangular lattice) en- 
ters our calculation via the dispersion of both c and / 
particles. Note that conventional magnetic order is ex- 
cluded from the outset due to the choice of the decoupling 
scheme. The lattice degrees of freedom are assumed to 
be 3d in all cases, but this does not explicitly enter the 
mean-field theory (6). The calculations are intended for 
T = 0, but are performed at a small finite T for conver- 
gence reasons. 

We obtain phase diagrams as function of external pres- 
sure and inter- moment exchange, with the bulk modulus 
B and the conduction band filling Uc as system parame- 
ters (nj = 1 is fixed in the Kondo limit). Most qualita- 
tive features are independent of the spatial dimensional- 
ity and underlying lattice geometry, therefore we discuss 
the 2d square lattice results unless otherwise noted. 

Sample phase diagrams for diflFerent values of the bulk 
modulus are shown in Fig. 2. A sufficiently soft lat- 
tice with a correspondingly small bulk modulus makes a 
first-order transition at the Lifshitz line preferable. Be- 
low a critical bulk modulus ^ this first-order lattice 
transition has no zero-temperature endpoint (Fig. 2a). 
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FIG. 3: As in Fig. 2, but for electrons on a 3d cubic lattice, 
for B = 0.005. 

Instead, it continues to Jh = which is the Kondo 
volume-collapse transition considered in Ref. 15. By in- 
creasing the bulk modulus above B*, the first-order tran- 
sition shows a quantum critical endpoint which is quite 
generically located on the Lifshitz line within numerical 
accuracy (Fig. 2b). Here, for larger values of the inter- 
moment exchange Jh^ the first-order transition coincides 
with the transition line to the fractionalized Fermi liquid. 
This type of transition, i.e., a direct first-order transition 
from a heavy Fermi liquid to a non-Kondo state domi- 
nated by inter-moment exchange (which is FL* here, but 
will be antiferromagnetic if frustration is reduced) also 
appears in other theories, but is usually not observed 
experimentally. By further increasing the bulk modu- 
lus, the first-order volume transition becomes weaker. In 
particular, the coincidence with the transition to the FL* 
phase is shifted to larger values of the parameter Jh (Fig. 
2c). Beyond a critical value of the bulk modulus, first- 
order transitions do no longer exist. Fig. 2d. This behav- 
ior turns out to be well described by a Landau theory, to 
be described in Sec. Ill below. 

As our calculations are performed at small finite T, 
care is required if the characteristic scales of the problem 
become small. If Jh or Tk become of order T, the tran- 
sition lines in the phase diagrams start to deviate from 
the zero-temperature solution. From numerical results at 
lower temperatures we obtained the asymptotic behavior 
of the transition lines which is indicated by the dashed 
lines in Fig. 2. 

The large- results for 3d electrons, moving on a cubic 
lattice, are qualitatively similar to those ior d = 2. An 
endpoint of the first-order lattice transition line is found 
(Fig. 3) - whether this endpoint can again be exactly 
located on the Lifshitz line without fine tuning is diffi- 
cult to decide, due to limited numerical resolution. (This 
question shall be investigated using the Landau theory of 
Sec. Ill below, with a negative answer.) Above a critical 
bulk modulus, the first-order volume collapse transition 
ceases to exist (not shown). Below another critical bulk 
modulus, the T = endpoint of the first-order transition 
disappears. As in = 2, the first-order line becomes a 
direct FLi-FL* transition for large values of the inter- 
moment exchange Jh- 

In our numerical results the mean-field bond ampli- 
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FIG. 4: Quasiparticle dispersions for different parameters 
along a constant- J/f cut of the phase diagram in Fig. 2, 
for parameters t = 1, ric = 0.8, Jh = 0.1, Jk = 1.5, 
7 = 0.05, B = 0.01 and T = 0. Upper panel: Bands along 
the — (7r,7r) (7r,7r) direction for pressures in the a) FLi, 
b) FL2 (close to the Lifshitz transition), and c) FL* phases 
{p = 0, 0.0085, 0.1, respectively). Lower panel: Compari- 
son of the low-energy part of the dispersions. Case a) has 
one heavy electron-like band intersecting the Fermi level. In- 
creasing pressure lowers the energy at the zone boundary and 
eventually causes the emergence of a hole-like sheet, case b). 
At higher pressures, the two bands evolve into conduction 
electron and spinon bands, which each intersect the Fermi 
energy once - this is the FL* phase, case c). 



tude xo is in all aforementioned cases negative, leading 
to a hole-like spinon dispersion (see also Fig. 4). At 
a Lifshitz transition, the lower band becomes sud- 
denly hole-like, promoting a hole-like pocket which splits 
from the Fermi surface. Obviously, the emerging hole- 
like pocket will contribute a power-law density of states 
at the Fermi energy, leading to a singular contribution to 
the thermodynamic potential - this will be discussed in 
more detail in the next section. 

A comment on the finite-temperature physics is in or- 
der: Both the Lifshitz and Kondo breakdown transitions 
continue to T > at mean- field level, however, it is clear 
that both are smeared into crossovers upon inclusion of 
fluctuation effects.^ The first-order lattice transition (for 
fixed p, Jh) ends in a finite-temperature critical end- 
point. These endpoints define a line in a generalized 
phase diagram which terminates in the T = endpoint 
described above. 

Let us summarize the physical picture obtained so far. 
Upon inclusion of lattice degrees of freedom, the Lifshitz 
transition between FLi and FL2 tends to become first 
order (because the system can gain energy by abruptly 
gapping one of the bands), while the Kondo breakdown 
transition remains of second order, except for soft lat- 
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tices where both transitions merge into a single first-order 
transition. The calculations further suggest an interest- 
ing coincidence of two critical phenomena in the case of 
2d electrons, namely the endpoint of the volume collapse 
is located precisely on the Lifshitz transition line. To 
understand the latter behavior in more detail, we now 
devise an effective theory for the two phenomena. 



III. LANDAU THEORY FOR LIFSHITZ AND 
LATTICE VOLUME TRANSITIONS 

Two prominent features of the heavy-fermion model of 
Sec. II were (i) Lifshitz transitions, where the Fermi sur- 
face topology changes with the appearance or disappear- 
ance of a Fermi sheet, and (ii) lattice transitions of the 
volume-collapse type. Note that Lifshitz transitions not 
only appear in the Kondo-breakdown scenario of Ref. 7, 
but also in scenarios for metamagnetic behavior inside 
the heavy Fermi liquid. These findings motivate to an- 
alyze an effective Landau-type theory for Lifshitz plus 
lattice physics - the mean-field analysis requires less nu- 
merical effort than that of the full fermionic large- A/" the- 
ory of Sec. II. The discussion of fluctuation effects will 
be postponed to Sec. IV. 



A. Landau functional 

The vantage point is a conventional Landau theory of 
an iso-structural volume-changing transition of the crys- 
tal lattice - this transition is of liquid-gas type, with 
a strain order parameter e, Eq. (2). Further, we assume 
that the system has a Lifshitz transition, where the chem- 
ical potential crosses a band edge. Importantly, the lat- 
tice parameter e couples linearly to the tuning parameter 
of the Lifshitz transition. (Microscopically, e changes e.g. 
hopping matrix elements, which in turn moves fermionic 
bands relative to the chemical potential.) 

The band edge crossing the Fermi level gives rise to 
a non- analytic piece in the Landau effective potential at 
zero temperature.^^ The fermion density of states in d 
dimensions close to the band bottom is: 



D(uj) (X 6(cj)|cj| 



(12) 



where Q{x) = (1) for x < (x > 0) is the step function. 
Then, the non-analytic part of the Landau free energy 
arising from the kinetic energy of the fermions is 



FF(e) = -/>:e(-e + eo)|e-eo| 



(13) 



where > is a coupling constant. For e < eo the band 
becomes occupied, where eo is the location of the Lifshitz 
transition. 

Choosing a reference volume Vq such that eo = 0, we 
can add a regular forth-order expansion in the order pa- 
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FIG. 5: (Color online) Samples of the mean- field free energy 
F{e) (14), for various values of the parameter v. The left 
figure corresponds to v > ^ and d = 2, the right figure to 
V < {) and d = 3. If m drops below mc2, a second minimum 
evolves at some e 7^ 0. If even m < md, this minimum 
becomes the absolute minimum (see also Tables I and II). 



rameter to obtain the Landau functional:^^ 



F{e) = Fo-Vop'^e- 
- ^e(-e)|e|^ 



m 



-e' + -e^ 



+1 



(14) 



where is a reference pressure. Microscopically, the 
mass m is related to the bulk modulus of the material, 
but will also depend on electronic properties (e.g. the 
ratio Jh/Tk in the Kondo lattice case). A cubic term e^ 
is symmetry- allowed and generically appears in the ex- 
pansion (14). (Note that in Ref. 15 the reference volume 
was chosen such that 'U = OatT = 0- this freedom is 
absent in the present case.) 

The physical situation corresponds to pressure tuning 
and is described by the free-enthalpy functional G(e) = 
F(e) +p]/be, which needs to be minimized to obtain the 
equilibrium value of the strain e. We shall discuss the 
arising T = phase diagram, upon variation of m and p, 
separately for the cases d = 2 and d = 3 in the following 
subsections. 

Before proceeding, we remind the reader that the ther- 
modynamic singularities at a Lifshitz transition are in 
general weak: The ground-state energy follows F ex 
-Q{h)\h\^ {F (X -e{h)\h\^/^) for electrons in d = 2 
(d = 3) as function of the tuning "field" h. In con- 
trast, at a usual second-order transition, the free en- 
ergy scales as F oc — \h\^^^ {F oc — 9(/i)|/ip) if h couples 
linearly (quadratically) to the order parameter, and we 
have assumed mean-field exponents (and neglected fluc- 
tuations). These considerations are relevant, because, for 
situations with a linear coupling between order parame- 
ter and pressure, it is known that lattice effects can ren- 
der the transition first order due to the compressibility 
diverging at a finite distance from criticality.^^ This effect 
originates from a diverging second derivative —d'^F/dh?^ 
and does not occur in the Lifshitz case. In our problem, 
the compressibility —d^G/dp^ will only diverge at the 
endpoint of the volume collapse transition. 
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d = 2 


0<\v 


< Vcl 


Vcl < \v\ < Vc2 


Vc2 < \v\ 




rrici 
mc2 


-— ef- 

9 u^y 
-— 6f- 

3 


-v) + 2k 
-v) + 2hi 


lir + 2«0(-«) 


lir + 2«e(- 


-v) 
-v) 



TABLE I: Location of critical points from Landau theory in 
= 2, in dependence of the cubic coefficient v. The critical 
values are Vd = V&ok, and Vc2 = Sy/uK,. If rrici < mc2, 
a critical endpoint mc2 away from the Lifshitz transition 
emerges. This is not possible if < ^; < t'ci since always 
rrici = mc2 = 2k, then. See text for details. 



B. Mean-field phase diagram in d = 2 

By construction, the Lifshitz transition in Eq. (14) oc- 
curs from e > to e < 0; a lattice transition may occur 
elsewhere on the e axis. For d = 2, the fermionic con- 
tribution in Eq. (14) is quadratic in e, i.e., renormalizes 
the mass m for e < 0. In the following, we analyze local 
minima of the free enthalpy G(e). We will obtain critical 
values mci,2 of the mass m: At m = rrici the Lifshitz 
transition (as function of p) changes from continuous to 
discontinuous, and at m = mc2 the lattice transition (as 
function of p) has its critical endpoint. 

A continuous Lifshitz transition is not possible if G(e) 
is not minimal at e = at any pressure. This is certainly 
the case if m — 2/^ < 0. If m — 2tz > 0^ the position of 
minima depends crucially on sign and magnitude of the 
cubic coefficient as depicted in Fig. 5. We distinguish 
these cases in successive order. 

(i) V > 0. Obviously G(e > 0) can acquire a second 
minimum if v becomes large such that 



dG(e) 
de 



p — p + me 



ve^ + ue^ 







(15) 



is possible if e > 0. The existence of two minima of G(e) 
is a necessary condition for a discontinuous transition 
to occur. The discriminant of the cubic equation (15) 
excludes a second minimum if v'^ < ^um. Obviously the 
first-order transition line terminates then at the critical 
endpoint = 2/^ as long as < ^' < 3^/uK:. In the 
opposite case, v > Sy^im, a second minimum of G(e > 0) 
develops for pressures that lead to three real solutions of 
Eq. (15). If V > |a/S^, in this second minimum G(e) 
becomes smaller than Fq at p = p"^ (Fig. 5, left). The 
Lifshitz transition changes from first to second order at 
^ci = giT' "^^^ first-order lattice transition line finally 
terminates at the critical endpoint mc2 = 3 17 where the 
second minimum vanishes (Fig. 6). Note that rrici = mc2 
if V < Vcl = V^UK, which leads to a critical endpoint 
located at the Lifshitz line (Fig. 6b). 

(ii) V < 0. Our analysis of positive coefficients v does 
not need to be modified if we repeat all arguments for an 
order parameter e < 0. In this case the Lifshitz transi- 
tion always becomes continuous at and 
the first-order lattice transition line terminates at the 



critical endpoint mc2 = I17 + 2/^. The phase diagram 
corresponds then always to Fig. 6a. 

(iii) V = 0. In this trivial case the critical endpoint 
rrici = mc2 = 2/^ is situated on the Lifshitz line. 

We summarize all different cases in two spatial dimen- 
sions in Table I. Qualitatively, all these cases divide 
into two different types of phase diagrams as depicted 
in Fig. 6. 

If the critical endpoint is on the Lifshitz line, then the 
effective potential at this special point is oc |ep (oc |ep) 
for positive (negative) e. This implies that e oc (p— p*)^/^ 
{e (X p* — p) for p > p* (p < p*) in the vicinity of the 
endpoint, i.e., for m = mci,2. The first-order line near 
this point follows |m — mci,2| oc |p — and along this 
line the jump in e is Ae oc |m — mci,2|. This shows that 
the critical exponents f3 and S do not take their mean- field 
values, due to the non- analytic behavior of the fermionic 
part of the Landau potential. 



C. Mean-field phase diagram in d = 3 

In three dimensions, the singular contribution to the 
free energy enters with a fractional power which is not 
contained in the series expansion of the regular part. 
In addition to the critical parameters rrici and mc2, 
it is necessary to introduce further critical parameters 
m'^^{u^ /^, v) and m'^2{^^ ^5 ^) which mark the onset of lo- 
cal minima of G(e) at e < (Fig. 5b). If pressure is used 
as control parameter, the relation of m to the critical val- 
ues mci, mc2, '^ci ^c2 determines the type of phase 
transition. If m < 0, generically a discontinuous Lifshitz 
transition occurs at a critical pressure p > p^ . Further- 
more, if -i; < we obtain ruci = mc2 = 0. By varying the 
parameter k and thereby tuning m^^ and entries 
of Table II can be realized. 

From Table II we can read off the phase diagrams 
possible in d = 3. Since always m'^i < zero- 
temperature endpoint of the lattice transition is never 
located on the Lifshitz line (except with further fine- 
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FIG. 6: Two representative mean-field phase diagrams in two 
dimensions, derived from the Landau theory. Thick/thin lines 
represent first/second order transitions. If the parameters 
rrici and mc2 coincide, the quantum critical endpoint of the 
lattice transition is located on the Lifshitz transition line (Fig. 
6b). In the other case, rrici 7^ mc2^ the first-order lattice 
transition line turns away from the Lifshitz line and ends at 
a critical endpoint (Fig. 6a). 
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d = S 


< m < rrici 


rrici < m < mc2 


mc2 < m 


< m < m'ci 


k LS 


1^' k LS " 


1'' & LS 


m'ci < m < m'c2 




-j^st b 




m'c2 < TTl 


& LS 




2nd 



^an additional 1^^ order transition can occur 
^two 1^^ order transitions occur 

TABLE IL Type of phase transitions from Landau theory 
in d = 3, for m > and v > 0. We denote a first- order 
transition which is also a Lifshitz transition by the entry 1^^ 
& LS. If m > and v < 0, we know from the 2d case that 
TTici — mc2 — and only the third column can be realized. 
The latter case leads again to the phase diagram of Fig. 6b. 
In the opposite case > we always obtain mc2 > tUci > 0. 
From the above table we conclude then that in addition to 
Fig. 6b, also the phase diagrams of Fig. 7 can be realized. 
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FIG. 7: Additional cases for mean- field phase diagrams that 
occur in d — 3. Both situations the first-order lattice tran- 
sition line bifurcates at a point which a) does [b) does not] 
coincide with the point where the Lifshitz transition changes 
from first to second order. 



tuning) . The behavior along the Lifshitz line can be again 
of the type shown in Fig. 6a, or there may be additional 
bifurcations of the first-order lattice transition (Fig. 7). 



D. Conclusions from Landau theory 

It has been shown that, for a range of mass parame- 
ters m, the Lifshitz transition is rendered discontinuous 
due to a coupling to the lattice, i.e., a first-order volume 
collapse coincides with the Lifshitz transition. The vol- 
ume collapse displays a quantum critical endpoint in the 
m-p plane. As a specific feature of the d = 2 theory, the 
volume-collapse endpoint is located on the Lifshitz tran- 
sition line in a large parameter regime. (In contrast, in 
d = 3 the Lifshitz and volume-collapse transitions sepa- 
rate before the volume- collapse endpoint is reached.) 



uous mean-field transitions remain of first order up to a 
critical temperature where a finite-temperature endpoint 
of liquid-gas type occurs.) 

The second-order Kondo breakdown transition has 
been investigated in Refs. 7 and 17 - this can lead to 
interesting non- Fermi liquid behavior, but we shall not 
repeat the arguments here. The Lifshitz transition is in 
the universality class of the dilute Fermi gas, with finite- 
temperature properties discussed e.g. in Ref. 34. 

It remains to consider the endpoint of the volume- 
collapse transition. This has a remarkable feature: 
The critical endpoint of an iso-structural volume-collapse 
transition has no critical fluctuations at finite wavevec- 
tors. To see this, first note that the strain field e is re- 
lated to longitudinal phonons by e V^; • 4> where is the 
displacement field. In continuum elasticity theory, it is 
well known that, even for a vanishing bulk modulus, a fi- 
nite shear modulus renders all phonon velocities finite. 
A detailed analysis by Cowley^^ shows that the same 
is true in crystals: All homogeneous strain fluctuations 
remain non-critical at the iso-structural volume- collapse 
transition. (Note that the vanishing of a phonon veloc- 
ity would imply a structural transition with symmetry 
breaking.) In other words, at the endpoint under con- 
sideration, only the zero-momentum mode becomes crit- 
ical, and Landau theory is exact for the lattice degrees of 
freedom. (Strictly speaking, phonon modes with wave- 
lengths of order sample size become critical.) 

For the interesting case of coinciding volume- collapse 
endpoint and Lifshitz transition, we can briefly discuss 
the thermodynamics. To this end, temperature correc- 
tions to the mean-field parameters and m of Eq. (14) 
have to be taken into account. Due to the underlying 
Fermi-liquid physics (of the additional bands crossing the 
Fermi level), these are quadratic in temperature, e.g., 
p* = p*(T = 0) + Ci^^- Furthermore, the fermionic con- 
tribution to Eq. (14) has to be replaced by the free energy 
of the Fermi gas. Near the Lifshitz transition this can be 
written as^^ 



FF(e,r) = -A^r2+i$(e/r) 



(16) 



which replaces Eq. (13), and ^{x) is a scaling function. 
Minimizing the free enthalpy G(e) in d = 2 now shows 
that e oc T upon raising temperature at the endpoint, i.e., 
the behavior is dominated by the fermionic part of the 
free energy. Consequently, the compressibility diverges 
as oc T~^, and the specific heat follows Cp oc T. 



DISCUSSION 



IV. BEYOND MEAN-FIELD THEORY 

So far, our analysis was restricted to mean- field the- 
ory derived either from a large- A/" or a Landau approach. 
In this section, we briefly comment on fluctuation effects 
near the continuous quantum phase transitions of our 
model system. (It is reasonable to assume that discontin- 



Our study of coupled Kondo-lattice and volume- 
fluctuation physics was motivated by antiferromagnetic 
quantum critical points in heavy-fermion metals where 
the proposal of a breakdown of the Kondo effect at crit- 
icality has spurred intense interest. Within the scenario 
of Ref. 7, where a transition between FL and FL* occurs 
via the vanishing of the quasiparticle weight on one of the 
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two sheets of the FL Fermi surface, we have investigated 
whether first-order volume transitions may spoil quan- 
tum criticality. These calculations merge former investi- 
gations on Kondo breakdown and Kondo volume collapse 
transitions. 

Our results show: (i) The Kondo breakdown transition 
remains intact as a second-order transition (except for 
very soft lattices), but the Lifshitz transition (i.e. a topo- 
logical splitting of the Fermi surface) which inevitably 
appears near the FL-FL* transition in the scenario of 
Ref. 7, tends to become first-order, i.e., occurs concom- 
mitantly with a Kondo volume collapse, (ii) Applied to 
the heavy- fermion metals CeCue-^cAua; and YbRh2Si2, 
this suggests to experimentally search for (possibly weak) 
first-order transitions at low temperature near the un- 
conventional candidate critical points. If those are not 
found, then alternative scenarios of Kondo breakdown 
should be investigated in more detail^' ^ - those may not 
necessarily imply the simultaneous existence of two Fermi 
sheets near the Kondo breakdown transition. Partic- 
ularly appealing seems the idea that, at the transition 
point, quasiparticles on both the large Fermi surface of 
FL and the small Fermi surface of FL* (or of the ordered 



state) become singular, i.e., that two hot Fermi surface 
sheets coexist, with a "super-large" Fermi surface of vol- 
ume 1 -h 2nc.^^ (As a result, the Lifshitz transition is 
removed from the agenda.) However, a mean- field the- 
ory for this scenario is not known, (iii) As an aside, we 
found that the zero-temperature endpoint coincides with 
the second-order Lifshitz transition for a large parameter 
regime in the case of 2d electrons. This coincidence of 
two critical phenomena, i.e. "non-Landau behavior", is 
rooted in the non-analytic free energy near the topologi- 
cal Lifshitz transition. 
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